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1 Introduction. 

The non-relativistic quantum models with position-dependent mass were used many 
years ago for the effective description of different models which are more complicated than 
the standard "academic" one-particle Schrodinger equation (as examples, see [1]). In this 
sense, position-dependent mass is frequently dubbed as effective mass (EM). During the last 
years this sort of models became again very fashionable in the literature, mainly due to the 
growing interest in such complex problems as theoretical descriptions of nanodevices, motion 
in curved spaces, and models with pseudo-Hermitian Hamiltonians (for illustration, see [2]). 
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Different approaches to this branch of modern Quantum Mechanics were explored, and 
super symmetric (SUSY) ones are among the most promising. In the framework of SUSY 
method, one may try to find functional dependence of effective mass and potential such 
that the corresponding Schrodinger equation is solvable or quasi-exactly (partially) solvable. 
These solvable models may be considered both as a base for perturbation expansions and as 
a laboratory for the study of qualitative properties. The SUSY method itself in Quantum 
Mechanics includes many different variants [3], and the SUSY intertwining relations seem to 
be the most promising among them. 

In most general form, SUSY intertwining relations between a pair of partner Hamiltonians 
if i j2 are: 



where the mutually Hermitian conjugate intertwining operators Q ± are called supercharges. 
These relations lead to the isospectrality of Hamiltonians Hi >2 up to possible zero modes of 
supercharges. This means that (again, up to zero modes of the energy spectra of Hi 
and H 2 coincide. Their bound state eigenfunctions are related (up to normalization factors) 
by the supercharges: 



H^\x)=E n ^\x); < = 1,2; n = 0,1,2,...; ¥^ = Q~¥^; = Q+< 2 >. (3) 



If either Q + or Q~ have some zero modes, and they coincide with the wave functions either 
of H 2 or of Hi, these wave functions are annihilated according to (3) and have no analogous 
states in the spectra of the partner Hamiltonian. Thus, from the intertwining relations one 
can find a pair of almost isospectral Hamiltonians, which sometimes can be solved (partially 
or exactly). By solving for the Hamiltonian we mean finding analytically its energy spectrum 
and corresponding wave functions. 

We stress that the intertwining relations approach described above is very general, pro- 
viding connections between pairs of spectral problems Hi ,H 2 . It depends neither on specific 



HiQ + = Q + H 2 ; 



Q~H X = H 2 Q 



(1) 
(2) 
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nature of spectral problems of the operators Hi, H 2 , nor on the specific form of intertwining 
operators Q ± [4], [5], [6], [7]. 

Each Hamiltonian Hi included in the intertwining relations (3) has at least one symmetry 
operator Ri = Q + Q~ , R 2 = Q~Q + , which commutes with it: 

[H i ,R i ] = 0; i = l,2 (4) 

Sometimes, these symmetry operators Ri are expressed in terms of the Hamiltonian Hi itself, 
and thus give no new information about the system. But otherwise, they describe indeed 
the symmetry of the model [6] , [7] . 

A variety of different realizations of the intertwining relations method was built during the 
development of SUSY Quantum Mechanics. In the original one-dimensional SUSY Quantum 
Mechanics with supercharges given by first order differential operators, not more than one 
normalizable zero mode of may exist, and the spectra of Hi, H 2 either coincide or differ 
by one bound state. The "symmetry operators" Ri coincide with the Hamiltonians in this 
case due to their factorization property: H x = Q + Q~; H 2 = Q~Q + . 

In multi-dimensional SUSY Quantum Mechanics developed in [8] the supercharges of first 
order in derivatives are used also, but a set of d related intertwining relations between two 
scalar and (d — 1) matrix Hamiltonians (d is dimensionality of space) must be considered 
simultaneously^. A more interesting realization of intertwining relations was constructed for 
one- and two-dimensional systems by using supercharges of second [4], [6], [7] (and higher [5]) 
orders in derivatives (Higher order or Nonlinear SUSY Quantum Mechanics). In particular, 
for these systems the number of zero modes of Q^ 1 may be more than one (their number 
and properties are under the control). In the case of d — 2, many integrable systems with 
nontrivial symmetry operators Ri were built [6] . Most of these systems are not amenable to 
separation of variables. 

The models with effective mass were also incorporated in SUSY intertwining relations 

approach [10], but until now, all these papers concerned one-dimensional problems. To our 

d Alternatively, intertwining of first order between multi-dimensional scalar Hamiltonians leads [9] to the 
systems amenable to separation of variables. 
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knowledge, the only work on higher- dimensional case was the paper of C.Quesne [11]. In this 
paper, the general system of equations equivalent to SUSY intertwining with scalar super- 
charges of first order in derivatives was illustrated by a particular two-dimensional model, 
where the effective mass depends on one of Cartesian coordinates only: M(x) = M(x\). In 
this case, the EM Schrodinger equation certainly allows the separation of variables in Carte- 
sian coordinates, and the initial problem reduces to a pair of one-dimensional Schrodinger 
equations. It would be interesting to extend this approach to arbitrary form of EM function 
M(x) and to higher orders of supercharges. 

The organization of the present paper is the following. In Section 2 the general intertwin- 
ing relations of first order for two-dimensional EM Hamiltonians are fully solved in terms 
of four arbitrary functions. The mass functions, potentials and coefficients in the intertwin- 
ing operators (supercharges) are built explicitly in terms of these arbitrary functions. We 
show that for the most general solution both intertwined Schrodinger Hamiltonians allow 
separation of variables in suitable coordinates 2/1,2/2- Section 3 concerns SUSY intertwining 
of second order. After general formulas are derived, several models with specific forms of 
supercharges are analyzed. The general solution for a particular structure of supercharge 
(the case of elliptic metric) is obtained. In Conclusions the separation of variables as well as 
the symmetry properties of the system (its integrability) and irreducibility of supercharge 
are discussed. 

2 General solution of first order intertwining relations. 

The two-dimensional Schrodinger operator with effective mass m(x) = m M(x) and 
h = 2m = 1 can be written in explicitly Hermitian form [12] as e 

H = ~[M a (x)diM p (x)diAP(x) + M' y (x)d i M l3 (x)d i M a (x)} + V(x), (5) 
2 

where x = (xi,x 2 ), <9« = d/dxi, the effective mass function M(x) is dimensionless function 
of coordinates, the real function V(x) is the potential, and a, /5, 7 are constant parameters 

e Everywhere below the summation over repeated indices is implied. 
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satisfying the physically motivated restriction a + (3 + 7 = —1. The specificparticular values 
for a, (3, 7 depend on the specific physical model which is described by the EM Schrodinger 
equation. 

It is convenient to rewrite (5) in the form: 

H = -d i jJ^a i + U(a!), (6) 
which unifies different choices of 01, (3, 7 in the additional terms of effective potential U(x) : 

Let us suppose that a pair of partner two-dimensional Hamiltonians Hi, H 2 of the form 
(6), both with effective mass M(x), and potentials Ui(x) , U 2 (x), satisfy the intertwining 
relations (1), (2) with most general first order supercharges: 

Q + = (Q-y = qi(x)d l + q(x); ^ = 1,2. (8) 

The explicit form of the intertwining relations (1), (2) seems to be too complicated in the 
two-dimensional case to be solved analytically. Only choosing some ansatzes for M(x) and 
V(x) one can simplify the problem to make it solvable: this was done in paper [11]. 

In order to find the general solution of intertwining relations one must make their form 
more tractable. First of all, it is useful to eliminate the first derivatives from (6) a suitable 
similarity transformation: 

h,2 = e^^E-^+^l); (9) 



A ( 2)M(f) 3[d k M(x)Ud k M(x) 
^ - U ^~^PW + V 2M3(f) " ( 1Q ) 

with some real function (j>(x), and 

q k (x)(d k M(x) 

q± = e -*WQ± e *(S) = qi (Z) dl + p (5); = + 2 V M(f) y . (11) 
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The function <p(x) can be determined from the condition of absence of first derivatives in 

d k M(x) , | 



= 2M(x) ' ^ = 2 lnM ^"> + Const - 



(12) 



This similarity transformation certainly is not unitary, and therefore, hi t2 are not Hermitian, 
and q + ^ {q~Y ■ Nevertheless, these operators satisfy the intertwining relation equivalent to 

(1): 

hq + = q + h 2 . (13) 

Equivalence means that solutions of (1) and of (13) are in one-to-one correspondence. Fur- 
themore, solution of (2) will be automatically obtained, being the Hermitian conjugate of 

(!)• 

Thus, we will study relation (13), whose solution seems to be technically much simpler 
than (1). It can be rewritten as a system of six nonlinear differential equations for the 
functions M, v i j2 , q iy p : 



M(x) 



d k qi(x) + (diq k (x) 



+ S ik q j (x)id j M(x)) =0; 



2M(x)v(x)q t (x) - 2(d iP {x)j - (A^f)) = 0; 



M(x) 



qi(x) diV 2 {x) ) - 2v(x)p(x) 



+ (A( 2 V(f)) = 0, 



(14) 
(15) 
(16) 



where all space indices take two values 1,2, and the function v(x) is defined as: 

2v(x) = Vi(x) — v 2 (x). 



(17) 



The general solution of the three Eqs.(14) can be derived as follows. After subtraction 
of equations for i = k, one obtains that d±qi = d 2 q 2 , i-e. q\ = d\F(x)] q 2 = d 2 F(x). Then 
equation for i ^ k gives that the real function F(x) satisfies A^F(x) = 0. Therefore, the 
coefficients qi,q 2 are combinations of analytical functions f(z) and f*(z*) of z = xi + ix 2 
and z* = x 1 — ix 2 , correspondingly: 



Ql (x) = i(f(z) - f*(z*)); q 2 {x) = f(z) + /*(**)■ 



18) 
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After that, the last equation of (14) takes the form: 

(fd z - fd z *)M{x) = -(/' - r')M(x). (19) 

Its general solution can be found by the substitution M _1 (x) = j(f(z)f*(z*))Q(x), because 
Q(x) must satisfy the first order equation: 

(fd z - rd z *){l(x) = 0. (20) 

Thus, the general solution of (14) consists of Eqs.(18) and: 

M -^-\^MIm + Ij^)- (21) 

where Q is an arbitrary real function of its argument. 

The last terms in both equations in (15) disappear due to (18). Then Eq.(15) takes the 
form: 

(d z *p(xf) = iM(x)v(x)f(z); (d z p(x)j = -iM(x)v(x)f*(z*). (22) 

Therefore, p(x) is expressed in terms of an arbitrary real function P of one specific real 
variable: 

^ 4 K !( /t|-/7^)»> (23) 

and the function v(x) has a compact form in terms of Q and P : 

* w - * - * - m ~ ! i^MI m + 1 7$y) • (24) 

where P' means the derivative in respect to its argument. 

It is evident now that, instead of the Cartesian arguments X\ ) X2 ) the alternative real 
space coordinates yi,y2 are much more convenient: 

f dz f dz* ( f dz f dz* \ , , 



/(*) y /*(**)' ~- u /(z) y /*(**)/" 

Finally, taking into account the previous results, the remaining equation (16) admits the 
general solution: 

v 2 (y) = fl( yi ) (p% 2 ) + P'(y 2 )) + r( yi ), (26) 



where a new arbitrary real function T is introduced. 

Summarizing the obtained results: we derived the general solution of intertwining rela- 
tions (13), represented in (18), (21), (23), (24) and (26) in terms of an arbitrary analytical 
function f(z) and three arbitrary real functions Q(yi), P(y 2 ),T(yi), restricted by physical 
reasons only: for example, Q must be positive. The use of variables y±,y 2 allows to rewrite 
the operators h± >2 from (9) in a rather compact form: 

hi,2 = -^(yi)Af +%/i)(p% 2 ) TPW) +T( yi ); Af ee d m d yi . (27) 
The intertwining operator q + from (18), (23) is also simplified essentially by use of yi,y 2 '■ 



q + = A(-d y2 +P(y 2 )), 



(28) 



i.e. it depends on y 2 only, having a form typical for one-dimensional SUSY Quantum Me- 
chanics with constant mass. The operators h\ )2 with "mass" 4f2 _1 (?/i) are intertwined (only 
in the variable y 2 ) by the operators q + , therefore their eigenvalues and eigenfunctions are 
interrelated. Performing the inverse similarity transformation to the physical Hamiltonians 
Hip, one can link their spectral properties by Q^ 1 as well. 

The variables y\ , y 2 are useful not only to express the most general form of intertwined 
Hamiltonians #1,2, but also to formulate the algorithm of explicit separation of variables 
in the spectral equations for operators hi >2 (analogues of physical stationary Schrodinger 
equations for Hamiltonians H 12 ) : 



-7^iK 2) + n( Vl ) (P 2 {y 2 ) =F P'(y 2 ) ) + r( yi ) 



(1,2) 



(29) 



4~~ Wi/ y 

Indeed, dividing (29) by |il(yi) one may search for its solutions, skipping for simplicity 
upper indices 1,2, as products 

i>n = r) n {y 1 )p n (y 2 ) (30) 

of solutions of two one-dimensional problems: 

1 



Pn{V2) = e n p n (y 2 ); 



r(yi) -E n )n- 1 (y 1 ) 



Vn(yi) = -e n r)n(yi)- 



(31) 
(32) 



Notice that the "potential term" in (32) depends explicitly on the "spectral parameter" 
of (31). Selecting the solutions rj n (yi), p n (y 2 ) of (32), (31), one has to remember that the 
conditions of normalizability may be changed by similarity transformation^: 

V n = eVn = eV(yiK(y 2 )- (33) 

Thus, solving two-dimensional Schrodinger equations with Hamiltonians Hi 2 by separation 
of variables for hi j2 , one has to find discrete values E n , such that the functions (33), in turn, 
constructed from solutions of two one-dimensional problems (32), (31), are normalizable. 

It is clear that the construction above starts from a choice of the analytical function 
f(z). In particular, it is instructive to consider the simplest case f(z) = f*(z*) = Const, 
where, because of (25), y± — x\, y 2 — x 2 up to a change of scale. In the same time, the 
mass M(x) depends on x\ only, due to the general solution (21). Thus, in this case the two- 
dimensional operators h^ 2 in (13) are intertwined by q + only in one variable x 2 , on which 
the mass does not depend M(x) = M(xi). The two-dimensional problem (1) becomes fully 
one-dimensional. The analogous result (up to replacing x\ <-> x 2 ) would be obtained for the 
choice f(z) = —f*(z*) = Const, where M(x) = M(x 2 ). 

Just the models with mass M(x), depending only on one variable X\, were considered in 
paper [11]. It is possible to derive the most general condition for functions f,f*,Q above, 
such that the mass M(x) does not depend on x 2 . From Eq.(21), the equation (d 2 M(x)^J = 
gives: 

(^))'= W %?_7gf W - (34) 
The r.h.s. has to be independent on y 2 as well, leading to the functional-differential equation: 

ff - rr - wf - rr + tr + (n 2 = o (35) 

with two possible solutions only: 

f(z) = az 2 + bz + c; (36) 
f( z ) = ae Xz + (3e~ Xz + 7. (37) 



■*It is clear from (12) that normalizability is not changed if the mass function M(x) is everywhere positive 
finite function. 
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Depending on the choice of constants, many different functions M(x\) are possible. As 
example: 

f(z) = az 2 ; yit2 = -^; n(y 1 ) = ay 2 1 ; M(x) = \ (38) 

QjZZ (XCC-^ 

Analogously (although more cumbersome), the choice f(z) = (az + b) 2 leads to one of the 
models of [11] with M(x) = 7/(0x1 + b) 2 . 

One might be interested in models with variables yi, y 2 , coinciding with polar coordinates 
p, 9. This condition leads to f(z) = bz, with b - an arbitrary real (or pure imaginary) 
constant. There is no need to satisfy (34) in this case: f2 remains an arbitrary function of 
yi, defining the form (21) of mass function. In particular, for real value of b the mass M(x) 
is a function of p only: M _1 (:r) = ^-p 2 fi(| hip). In full analogy to the case of Cartesian 
coordinates above, the most general condition, when mass depends only on one of polar 
coordinate, is wider than the simplest form f(z) = bz. 

3 Second order intertwining. 

As it was shown in the previous Section, the first order intertwining of the two- 
dimensional operators hi^ leads actually to separation of variables in terms of yi,y2- And 
furthermore, after separation one obtains the intertwining problem (13) with a constant mass. 
This rather trivial result is not very surprising in the light of the study of two-dimensional 
SUSY Quantum Mechanics (with constant mass), where nontrivial achievements were ob- 
tained following two directions. The first one [8] dealt with two-component (vector) first 
order supercharge operators, which intertwined a given scalar Hamiltonian with matrix part- 
ner. Alternatively, if one doesn't wish to deal with matrix potentials, a second option was 
proposed [6] - intertwining by scalar supercharges q + of second order in derivatives. This 
latter approach leads us to use the second order operators in the problem with effective mass 
considered in the present paper. 

The most general intertwining operator of second order in (13) is: 

q + = gi j{x)didj + Ci(x)di + B(x), (39) 
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where gij, Ci, B are arbitrary real functions of x. Equating coefficients of each power in the 
derivatives in (13), one obtains a very complicated system of nonlinear differential equations 
for functions M(x),v li2 (x) and all coefficient functions in (39). In order to have a chance to 
solve this system, we make a simplifying ansatz for the metric in the supercharge: 

g ik (x) = Lj(x)5 ik . (40) 

From third order terms in (13) it follows that the mass and metric are related: 

M(x) =u'\x). (41) 



The second order terms give a system of three equations: 



M(x) (diCkffl + dkdig) 



2M 2 (x)v{x)u{x) - Cj{x) djM(x) 



(42) 



where v(x) was defined in (17). Two of these equations can be solved in terms of an arbitrary 
analytical function g(z) similarly to (18): 



d(f) = i(g(z) - g*(z*)); C 2 (x) = g{z) + g*(z*). 



(43) 



The solution of the remaining equation in (42): 



diCi(x) ) — v(x) 



2M(x) 



we postpone for later. 

The first order terms can be written as: 



+ C i (x)(d i M(x)^j =0. 



(44) 



di\B(x) + v 2 (x) = M(x)v(x)Ci(x). 



(45) 

Instead of the second equation in (45), it is more convenient to solve the special combination 
of both equations, obtained by multiplying by Ci and subtracting: 

C 2 (f)9! - d(f)9 2 ) (b(x) + v 2 (xf) = 4d tl (b(x) + v 2 (xf) = 0, (46) 
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where we used the variables ti^ in full analogy to y^ 2 , defined in (25), up to the replacement 
f(z) <-> g(z). Therefore: 

B(x)+v 2 (x)=4S(t 2 ), (47) 

with arbitrary function S(t 2 ). After that, we express d\ in the first equation of (45) in terms 
of d tl ,d t2 and use (47), deriving the relation between mass M(x) and v(x) : 

V ; g(z)g*(z*)v(x) V ; 

Now we can solve the postponed Eq.(44), rewriting it in terms of the variables t 1 ,t 2 
and using (43) and (48). After some transformations, it takes the form of the well known 
Bernoulli equation for the function v(t): 

d ^i) = \v\t) + v(t)^M. (49) 

Its general solution [13] (section 1.1.5) includes an integration constant, which becomes an 
arbitrary function A(ti) of ti in our case: 

The term without derivative operators in (13) reads: 

A {2) (b(x) + v 2 (x)) + M(x)d(x) (diV 2 (x)) - 2M(x)v(x)B(x) = 0, (51) 



and it can be transformed to the first order differential equation for v 2 (ti, t 2 ) with coefficients 
which were defined above: 



S'(h 



2i K S'(t 2 ) + S\t 2 )J 



(dtMt)) + 7 — vMi) ~ -7 r~ = 0. (52) 

(Aih) + S(t 2 )j (Aih) + S(t 2 )j 



The general solution for v 2 and its partner v± have the form: 

2 



Vl 



At) 



S(t 2 )+A(t!) 



(53) 
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where D{t\) is a new arbitrary function. Finally, the operators take the form: 

2 . ™ 2 



hl ' 2 5(t 2 )+A(t 1 ) A » ) + 5(fo) + .4(A i ) 



(54) 



allowing the separation of variables t 1 and t 2 in the corresponding Schrodinger equations, 
analogous to (29) of previous Section. But in contrast to Eqs.(31), (32), the spectral param- 
eter E n is present now in both one-dimensional equations, leading to essential difficulties in 
solving the spectral problem. 



4 Conclusions. 

As we discussed in the Introduction, an arbitrary interwining relation between phys- 
ical (Hermitian) Hamiltonians Hi 2 leads to the Hermitian symmetry operators R\ = 
Q + Q~, R 2 = Q~Q + ■ In the context of Section 2, these operators are obviously of second 
order in derivatives. Meanwhile, in Section 3, they are initially of fourth order. Since the 
symmetry operator is defined up to a function of Hamiltonian itself, one must study the 
possibility to reduce the order of the symmetry. Because from the very beginning we re- 
stricted ourselves to real coefficient functions in the supercharges Q ± , the terms of odd order 
in derivatives, being non-Hermitian, cannot appear in As for the term of fourth order in 
derivatives, due to relation (41) between metric of supercharge and mass function in ifj, it 
can be just expressed as Hf plus some terms of even (second and zero) orders. This means 
that for the metric (40) in the supercharge, systems with Hamiltonians Hi obey symmetries 
of second order in derivatives only. This fact is not surprising in the light of separation of 
variables in Schrodinger equations discussed in the very end of Section 3. 

An interesting question concerns also the reducibility of the considered second order 
supercharges, i.e. the possibility in principle to factorize the operator q + onto two multipliers 
of first order of the form (11). It is useful to compare the terms of highest (second) order in 
(39), (40) and in (11), (18). After a simple analysis (in particular, due to absence of mixed 
term <9i<9 2 in (39)), it becomes clear that such factorization is possible only in the case when 
both the functions fi(z), f2(z) in the first order multipliers and uj{x) in (40) are constants. 
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This situation corresponds to the case of constant mass M(x) and is not relevant for the 
present discussion. Thus, the supercharges q + of the form (39), (40) are irreducible. The 
same conclusion can be reached after comparing the general expressions for the intertwined 
Hamiltonians (27) and (54) for the first and second order intertwinings. The denominator 
^S(t<z) + Aftx^J, present in both kinetic and potential terms for second order intertwining, 
may appear in the expressions (27) for Hamiltonians, involved in two consecutive first order 
intertwinings, for the case of constant mass only. This argumentation corroborates the 
conclusion about irreducibility of second order q + from another point of view. 
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